Exact valence bond entanglement entropy and probability distribution 
in the XXX spin chain and the Potts model 
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By relating the ground state of Temperley-Lieb hamiltonians to partition functions of 2D statis- 
tical mechanics systems on a half plane, and using a boundary Coulomb gas forrnalism, we obtain 
in closed form the valence bond entanglement entropy as well as the valence bond probability dis- 
tribution in these ground states. We find in particular that for the XXX spin chain, the number 
A'c of valence bonds connecting a subsystem of size L to the outside goes, in the thermodynamic 
limit, as {Nc){^) = ^ InL, disproving a recent conjecture that this should be related with the von 
Neumann entropy, and thus equal to jyK2 ^^L. Our results generalize to the Q-state Potts model. 

PACS numbers: 03.67.-a 05.50+q 









I 

o 
o 



m 
en 



o 



X 



Introduction. Entanglement is a central concept in 
quantum information processing, as well as in the study 
of quantum phase transitions. One of the widely used 
entanglement measures is the von Neumann entangle- 
ment entropy S'vNj which quantifies entanglement of a 
pure quantum state in a bipartite system. To define S'vn 
precisely, let p = |^)('I'| be the density matrix of the sys- 
tem, where |^) is a pure quantum state. Given a com- 
plete set X of commuting observables, let X = AU B he 
a bipartition thereof. Then S'vn is defined as 



S'vn(^) == -TlA PA In PA 



(1) 



where pA = Trs p is the reduced density matrix with 
respect to A. One readily establishes that S'vn(^) = 
Sv^iB). In most applications, the subset A corresponds 
to a commuting set of observables characterizing only a 
part of the whole system, and so one may think of A 
denoting simply that subsystem. 

Critical ground states in ID are known to have en- 
tanglement entropy that diverges logarithmically in the 
subsystem size with a universal coefficient proportional 
to the central charge c of the associated conformal field 
theory. Let L = \A\ be the size of the subsystem, and 
N the size of the whole system, both measured in units 
of the lattice spacing, with I -^ L ^ N (as we shall 
invariably assume in what follows) . Then [l|, |2| 



SvN(A) = (c/3)lnL, 



(2) 



where by = we denote asymptotic behaviour. Away from 
the critical point, the entanglement entropy saturates to 
a finite value, which is related to the correlation length. 
The von Neumann entanglement entropy is not an easy 
quantity to calculate, analytically or numerically; nor is 
it easy to grasp intuitively. The plethora of algebraic 
and geometric reformulations of quantum hamiltonians 
in one and higher dimensions [3|, |j| suggests that there 
should be a more convenient, geometric way to define an 



entanglement entropy. In the case of systems admitting 
infinite randomness fixed points, in one [5|,|6|,|7[ as well as 
in higher [8|, \^ dimensions, the ground state \Q) can be 
represented as a single valence bond state, and S'vn coin- 
cides with the number of singlets that cross the boundary 
of A times the logarithm of the number of states per site. 
In an interesting recent paper [1(| (see also 111 for a 
related, independent work) it was suggested that, even 
when \fl) is not a single valence bond state but a super- 
position of such states, the average number of singlets 
A'c(ri) crossing the boundary of the subsystem (multi- 
plied e.g. by In 2 for spins 1/2) could still be used as 
a measure of the entanglement entropy with all the re- 
quired qualitative properties. Moreover, it was observed 
numerically in 10] that, up to statistical errors, this 
valence bond entanglement entropy Svb bad the same 
asymptotic behavior as Svn for the XXX quantum spin 
chain, namely ^ with c = 1. In other words, the obser- 
vation of idl was that 



{Nc)m 



3 log 2 



InL- 0.481 InL. 



(3) 



Apart from this, the valence bond basis has been actively 
studied recently, in particular from rigorous |12l . |l3[ and 
probabilistic [ij] points of view. 

Studying entanglement entropy from a geometrical 
point of view seems particularly appealing in view of the 
quantum dimer model of Rokhsar and Kivelson Il5| and 
its many recent generalizations. Indeed, Henley [16| has 
shown that for any classical statistical mechanics model 
equipped with a discrete state space and a dynamics sat- 
isfying detailed balance, there is a corresponding quan- 
tum Hamiltonian whose ground state \n) is precisely 
the classical partition function. Of particular interest 
are then statistical mechanics models in which the mi- 
croscopic degrees of freedom directly define the valence 
bonds. This is the case for a certain class of lattice mod- 
els of loops, to be studied below. 
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FIG. 1: Loop representation of the periodic TL algebra (with 
A^ = 8) on the infinite half cylinder C-. The basis state 
corresponding to the upper rim is •(())()•. 



We show in this Letter that the probabihty distribution 
of the number of singlets crossing the boundary can be 
exactly determined for the XXX spin chain as well as for 
the related Q-state Potts model hamiltonians. We find 
that (m is not quite correct: the exact leading asymptotic 
behaviour is in fact {Nc){Vl) = ^InL ~ 0.405 In L. AU 
other cumulants have similar closed form expressions. 
Entanglement and the TL algebra. The 2D classical 
Q-state Potts model can be defined for Q non integer 
through an algebraic reformulation where Q enters only 
as a parameter. For this, recall that the transfer matrix 
in the anisotropic limit gives rise to the hamiltonian [22| 



H 



N-l 



(4) 



Here the e^ are elements of an associative unital algebra 
called the Temperley-Lieb (TL) algebra, defined by 



(5) 



[e„e,] = Ofor |i-j|>2 



The ground state of H depends on the representation 
of the algebra. For our purpose, it is natural to use 
the loop model representation, where the generators act 
on the following non-orthogonal but linearly independent 
basis states. Each basis state corresponds to a pattern 
of N parentheses and dots, such as () • (())•. The paren- 
theses must obey the typographical rules for nesting, and 
the dots must not be inside any of the parentheses. These 
rules imply that the () pairs consist of one even and one 
odd site, and that dots are alternately on even and odd 
sites. We start by convention with an odd site on the left. 
Note that we have used an open chain for convenience, 
but that a periodic chain can be considered as well. This 
requires the introduction of an additional generator cn 
coupling the N^^ and first site, and in the graphical rep- 
resentation, parentheses can now be paired cyclically so 
)•(())•( is now a valid pattern. This is shown in Fig. [TJ 

One can interpret these states in terms of spin 1/2 
by associating with each pair of nearest parentheses () a 



Uqsl{2) singlet {y/Q = q + q~^) so a valence bond can be 
drawn between the two corresponding sites. (Note that 
the generator a is nothing but the operator that projects 
sites i and i + 1 onto the singlet.) For the dots, the state 
must be chosen such that the application of the projec- 
tion operator onto the UqSl(2) singlet for any two dots 
that are adjacent (when parentheses are ignored) annihi- 
lates the state. Thus those sites are "non-contractible" . 
With these definitions, it is clear that the TL algebra does 
not mix states with different numbers of non-contractible 
sites. For Q generic, the set of basis states with fixed 
number 2j of such sites provides an irreducible represen- 
tation of TL, of well-known dimension 



N 
N/2 + j 



N 
N/2 + j + l 



(6) 



where N/2 + j must be an integer. This dimension co- 
incides with the number of representations of spin j ap- 
pearing in the decomposition of the product of N spins 
1/2. This is no accident: it is well-known [17] that the 
uncrossed diagrams are linearly independent and form a 
basis of the spin j sector in the sl{2) case; the results 
extends trivially to the Uqsl{2) case with q generic. 

When ^/Q > the ground state is found in the sector 
with j = for N even (and j = 1/2 for N odd). Note that 
the valence-bond basis is not orthonormal. The simplest 
way to proceed is thus not to calculate matrix elements 
of the hamiltonian H in this basis {wi\H\wj) but rather 
to define a non-symmetric matrix hij by expressing the 
action of H on any state as a linear combination of states 



j 



(7) 



The matrix hij is unique due to the linear independence 
of the states. The eigenvalues and right eigenvectors of 
h give those of H. 

Since all entries hij are (strictly) positive, the Perron- 
Frobenius theorem implies that the ground state jfijex- 
pands on the basis words with positive coefficients [23| 



|f})^^A, 



A^ > 



(8) 



We define the number of valence bonds Nc connecting 
the subsystem to the outside as the number of unpaired 
parentheses in the subsystem. We are here interested in 
its mean value 



{Nc)m 






and more generally in the probability distribution 



^{Nc) 



Y^w^-u 



(9) 



(10) 



Below we establish the leading asymptotic behaviour of 
(Nc) (and the higher cumulants) in the scaling limit 1 ^ 



L -^ N. Note that the TL formulation shows relationship 
between the Potts hamiltonian when ^/Q — 2 cos -j^ , 
with k integer, and the interacting anyons (coming in 
k + 1 species) hamiltonian in [18|. The valence bond 
entanglement entropy can be defined for these models as 
well, and, in the sector of vanishing topological charge, 
coincides with the one we are studying. 
Mapping onto a boundary problem. The wave func- 
tion in the ground state of a certain hamiltonian with pe- 
riodic (free) boundary conditions [2J] can be obtained as 
the path integral of the equivalent Euclidian theory on a 
infinite half cylinder (annulus) , denoted C_ (or A^ ) . To 
translate this in statistical mechanics terms, note that if 
we consider the square lattice with axial (diagonal) direc- 
tion of propagation [cf. Fig. [T], the hamiltonian belongs 
to a family of commuting transfer matrices describing the 
Q-state Potts model with various degrees of anisotropy. 
The ground state of all these transfer matrices is given 
by |0). Let us chose for instance the particular case 
where the Potts model is isotropic, with coupling con- 
stant e^ = 1 + ^/Q. Now the ground state \il.) can be 
obtained by applying a large number of times the trans- 
fer matrix on an arbitrary initial state (corresponding to 
boundary conditions at the far end of C_ or A_ .) Clearly, 
by the mere definition of the transfer matrix, this means 
that the coefficients of the ground state jfi) on the words 
\w) are (up to a common proportionality factor) equal to 
the partition function of the 2D statistical system on C_ 
or j4_ with boundary conditions specified by jw). 

We must now study such partition functions. We move 
immediately to the limit N ~* oo. We then have a system 
in the half plane, which, in the geometrical description, 
corresponds to a gas of loops with fugacity y/Q in the 
bulk, with half loops ending up with open extremities on 
the boundary. To go to the continuum limit it is conve- 
nient to transform this loop model into a solid-on-solid 
model [19|. For this parametrize ^/Q = 2cos7reo with 
< eo < 1 [25( Give to all the loops an orientation, 
and introduce complex weights exp(±^^) for the left 
and right turns. Since on the square lattice the num- 
ber of left (iil) minus the number of right (nn) turns 
equals ±4, this gives closed loops the correct weight ^/Q. 
Meanwhile, loops ending on the boundary will get, with 
this construction, the weight -^/Qb = 2 cos ^^ since for 
them riL — nu = ±2 [20|. Although no such boundary 
weight appeared in the initial lattice model and partition 
function, we note that for the fully packed loop model we 
are interested in it does not matter: the number of open 
loops touching the boundary is just N/2 — j, a constant. 

Introducing this boundary loop weight allows complete 
mapping to the SOS model (or free six- vertex model). 
Now it is known that in the continuum limit, the dy- 
namics of the SOS height variables turns into the one 
of a free bosonic field [19| . In a renormalization scheme 
where loops carry a constant height step A(f> = ±7r, the 
propagator of the field evaluated at two points x, x' on 



the boundary reads, in the infinite size limit [2C 

($(x)$(a;'))N 



-—Inlx — x'\ 



fll) 



where g — 1 — eo. Here the subscript N indicates Neu- 
mann boundary conditions, corresponding to the pres- 
ence of loop extremities on the boundary. 

Let us now single out a segment of length L on this 
boundary and attempt to count the number of loops con- 
necting this segment to the rest of the boundary. To do 
this we insert a pair of vertex operators, one at each a ex- 
tremity of the segment, V = exp [i(±ei -I- eo/2)$]. These 
operators do not affect the loops encircling the whole in- 
terval L since they modify the weight of such loops from 

g±i7reo/2 ^o ^±i7reo/2 ^Ti^co ^ gT*'reo/2^ -^l^^g giving the 

same sum \/Qb- ^^^ ^^r loops connecting the inside to 
the outside, the weight is now w = 2cos7rei. The bound- 
ary dimension of the fields V is, using the propagator 

4e? - en , , 

/i=^ (12) 



so their two-point function decays as L~^^. We can then 
find the average number of loops separating two given 
points by taking appropriate derivatives, and setting ei — 
eo/2 in the end. This leads to our main result 

2cos(7reo/2) 



{Nc)m 



eo 



■InL 



(13) 



7r(l — eo) sin(7reo/2) 

For the XXX chain (eo = 0) this reads (iVc) = ^\nL^ 
0.405 InL, while for bond percolation (Q = 1 or eo = 1/3) 
we have (Nc) = — InL « 0.551 InL. The slope becomes 
1 exactly as eo -^ 1, or ^/Q -^ —2. We note that the 
result for the XXX case is close but definitely different 
from the one proposed in [lO|. 

It is amusing to observe that one can exactly interpret 
the valence bond as singlet contractions for an ordinary 
supergroup in the case Q = 1, by taking a lattice model 
where the fundamental three-dimensional representation 
of SU{2/1) and its conjugate alternate. The hamiltonian 
is again ^, but this time the Ci are projectors onto the 
singlet in 3 (g) 3. The effective central charge for the this 
spin chain is CoS = 1 + -T-[arccosh(3/2)l^ ~ 1.845, and 
extending the argument suggested in [Ifl] for the XXX 
case gives a slope of j^f^ ~ 0.559, even closer to the 
exact result (IT^ . 

Of course, by taking higher derivatives of the two-point 
function of the vertex operators one can access the higher 
moments of (jlOp . In fact, the two-point function itself 
is nothing but the characteristic function oi p{Nc), al- 
though carrying out the Fourier transform in general is 
somewhat cumbersome. We will content ourselves here 
by giving the first few cumulants Ck = (cfc/7r'^') InL, with, 
in the XXX case (top) and the Q = 1 case (bottom): 

_ 4 _ / 8/3 _ r 16/15 

"^'^ V^TT ""' " \ 2(2^V3 - 9) ""^'X 8{5ttV3-27) 

(14) 




FIG. 2: Comparison between exact and numerically deter- 
mined values of the slopes ci and C2, shown as functions of 
the parameter eo. 



-2, the 



together with the observation that, as \fQ 
probabihty distribution becomes Poissonian: 



lim P(7Vc 



-InL 






N, 



(15) 



Numerical calculations. We have computed the dis- 
tribution (jlO[) numericaUy by exactly diagonalizing the 
transfer matrix, for periodic chains of size up to iVmax = 
32. The cumulants Ck on InL of p{Nc) obey a very sim- 
ple finite size scaling (FSS) form, where InL has to be 
replaced by -^ In sin -^ ; this follows from standard for- 
mulas for two-point functions of our vertex operators V. 
Precise values of the slopes Ck can then be extracted from 
a careful analysis of the residual FSS effects. As shown in 
Fig. [2] they agree well with our analytical results, except 
for Q ^ 4, where we expect logarithmic FSS corrections. 
For (5 = 1, the combinatorial nature of \Q) implies 
that all A^ in ^ are integers. This allows to obtain 
p{Nc) exactly for finite L and TV < iVmax- Using this, we 
can in some cases conjecture p{Nc) for any value of N 
[2lj . In particular we have established that 



N/2-1 

(N,) ^ {N^ ^ L^)pk{N') H {N' 



(2n+l)2)"-^/2 



(16) 
where pk is a polynomial of degree k = ^{N + 4:){N — 2) 
in N'^. This exact FSS form allows to obtain for the slope 
ci — 0.5517 ± 0.0003, in very precise agreement with the 
value 0.551329 from (fT3l). 



Conclusions. As already argued in [lOj, S'vb is a mea- 
sure of entanglement which seems as good qualitatively as 
S'vN, and easier to obtain numerically. We have shown in 
this Letter that it is also possible to tackle it analytically 
in the ID case. The results are less easily expressed than 
for S'vN (which is proportional to c) . On the other hand, 
they fit considerably more naturally within the trans- 
fer matrix and Coulomb gas formalism. It remains to 



be seen what happens for other models, and whether in 
particular a c-theorem of sorts is obeyed for S'vb ■ 
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The scale of Q affects the sound velocity and is impor- 
tant when studying the scaling of gaps. But it does not 
matter when dealing with entanglement issues. 
Note that the scalar product of a state \w) with itself 

is equal to ^/Q ^ for all w, so there is no need to 

consider "normalized" basis words. 

The boundary conditions at large N are not expected to 

affect the leading behaviour of S'vb. 

When ^/Q < (i.e. i < eo < 1) the true ground state has 

spin j > 0. Nevertheless, we continue to let \i}} denote 

the j = ground state. Numerical studies then indicate 

that ^ holds even for ^/Q < 0, if A is large enough. 



